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Abstract. The two-fold aim of the paper is to unify and generalize on the one hand 
the double integrals of Beukers for £(2) and £(3), and of the second author for Euler's 
constant 7 and its alternating analog ln(4/7r), and on the other hand the infinite products 
of the first author for e, of the second author for ir, and of Ser for e 7 . We obtain new 
double integral and infinite product representations of many classical constants, as well 
as a generalization to Lerch's transcendent of Hadjicostas's double integral formula for 
the Riemann zeta function, and logarithmic series for the digamma and Euler beta func- 
tions. The main tools are analytic continuations of Lerch's function, including Hasse's 
series. We also use Ramanujan's polylogarithm formula for the sum of a particular series 
involving harmonic numbers, and his relations between certain dilogarithm values. 
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1. Introduction 

This paper is primarily about double integrals and infinite products related to the Lerch 
transcendent <&(z, s, u). 

Concerning double integrals, our aim is to unify and generalize the integrals over the 
unit square [0, l] 2 of Beukers jl] and Hadjicostas P, [TU] for values of the Riemann zeta 
function (Examples 13.81 and 14. lj) . and those of the second author [201, f° r Euler's 
constant 7 and its alternating analog ln(4/7r) (Examples 14.41 and I3.14j) . We do this in 
Theorems 13.11 and 14. 1[ using a classical integral analytic continuation of $ (Lemmas 
12.11 and 12. 2j) . As applications, we obtain new double integral representations of many 
constants, including ln2, hup, n 3 , lncr, T(3/4), C(5), and G/n (Examples I3.4[ ETKj 13.91 
13.131 13.191 13. 2H and I3.24|) . where tp is the golden ratio, a is one of Somos's quadratic 
recurrence constants, and G is Catalan's constant. (Single integrals for them follow by the 
change of variables in the proof of Theorem 13.11 ) Corollary 13.11 gives a double integral for 
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the polylogarithm function; examples use Ramanujan's relations between certain values 
of the dilogarithm (Example IH.fi)) . Corollaries 13.5) RTTH1 and l4~Tl evaluate certain types of 
double integrals in terms of the function lnT(a;) and its derivative, the digamma function; 
Examples 13. 14) 13. lfil and 14.41 for ln7r, ir/\^3 and 7 are interesting special cases. Theorem 
13 d evaluates a class of double integrals with the help of Ramanujan's summation of a series 
involving harmonic numbers; applications use classical values of the di- and trilogarithm 
(Example EM) - 

Regarding infinite products, our aim is extend those of the first author (see [231) f° r e 
( Example 15.12)) . of the second author [23J for n ( Example 15 .1)) . and of Ser [T7j (see also 
|2*T] . [23J) for e 7 ( Example 15.8)) . To do this we use two analytic continuations of $ by 
series, one due to Hasse [12] (Theorems 12.11 and 12.2)) . We find new infinite products for 
many constants, among them r(l/4)r(3/4) _1 , e G ^ w , n/e, e n , a, and the Glaisher-Kinkelin 
constant A ( Examples 15.4) 1531 15.fil 1531 15.101 and 15.11)) . The product for a has rational 
factors. 

Finally, we mention other results obtained: explicit formulas for the Bernoulli and Euler 
polynomials ( Examples 12.41 and 12.5)) . and logarithmic series for the digamma and Euler 
beta functions (Theorems 15.11 and 15.2)1 . The latter series leads to infinite products for 
the exponential function of the same shape as other products in the paper, but which we 
construct without using Lerch's function (Theorem 15.31 and Examples 15. 1 21 and 15. 13)) . 



2. The Lerch Transcendent 

The Lerch transcendent $ (see [21 section 1.11], |2~4") section 64:12]) is the analytic 
continuation of the series 

1 z z 2 

W H^s,u) = — + 7 , Ua + , +-" » 

u s [u + l) s [u + 2) s 

which converges for any real number u > if z and s are any complex numbers with 
either \z\ < 1, or \z\ = 1 and dt(s) > 1. Special cases include the analytic continuations 
of the Riemann zeta function 



00 

(2) c( S ) = E^ = $ ( 1 ' s ' 1 )' 

the Hurwitz zeta function 



k s 
k=i 



k=0 



= $(l,s,u), 

u + k) s 



the alternating zeta function (also known as Dirichlet's eta function 77 (s)) 

k s 



(3) C(s) = J2 [ -i^ = H-l,sA) 

k=i 

the Dirichlet beta function 
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the Legendre chi function 



x y 2k+l 



(5) Xs( Z ) = J2^-- = 2-^Us,- 



(2k 

k=0 v 



the polylogarithm 



oo u 



(6) Li n (z) = ^- = ^(z,n,l) 



k=l 



and the Lerch zeta function L(\, a, s) = $(exp(27rzA), s, a). 

From (P), Lerch's transcendent satisfies the following identities: 

(7) $(*,s,it + l) = -($(*, s, it)- — 

2 \ li s 

(8) $(z,s- l,it) = ( U + ^J ®{z,s,u), 

1 9$ 

(9) $(z,s + l,u) = — — (z,s,u). 

s aw 

Henceforth, we assume that u and v are positive real numbers. The following two 
lemmas are classical. (See J2|, [IS], (2IJ, and |2*5].) 

Lemma 2.1. 7/ either \z\ < 1 and 3?(s) > 0, or z = 1 and 3?(s) > 1, then 

do) $(z , s ,„ ) = _y o — t »-v t . 



Proof. Replacing 6 t",^* with ^ fc>0 ^ fc e ( u + fe )* ; and noting that 



(u-l)t 
e t —z 

oo 1 />oo 



11) / e-^Vd^Arf^-Y- 1 ^ J ' ( ' S) 



(u + k) s J Q ' ' (n + /c) s ' 
the lemma follows. □ 

Lemma 2.2. TTie series (0) extends to a function <3>(z, s, it) which is defined, and holo- 
morphic in z and s, for z G C — [l,oo) and all complex s, and which is given by the 
integral formula (JTUJ) i/3?(s) > 0. 

Proof. If z G C — [l,oo), the denominator of the integrand in (jl()j) does not vanish. It 
follows that for dt(s) > the right side of (JTUJ) defines a function which is holomorphic in 
z and s. The extension to other values of s follows inductively using the identity (JBJ). □ 

For z G C — [1, oo) and s = 0, —1, —2, . . . , there is a closed formula for the function 
s, it). In fact, summing the series $(z, 0, it) = 1 + z + + • • • and using (JBJ), we see 
inductively that the functions 

(12) $(2,0,«) = ^, 
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are rational, with rational coefficients, holomorphic in z except for a pole at z = 1. In 
particular, setting z — — 1 and m = — s, we define the mth Euler polynomial 

(14) E m {x)=2Q{-l,-m,x). 

On the other hand, if s and u are positive integers, then Q(z, s, u) is a transcendental 
function. For example, when s = u — 1 series (JTJ) is 1 + | + ^ + • • • , which evidently 
sums to 

(15) $0,M) = — ln(l-z) = -Lii(«). 

For s = 2, setting it = 1 in (jHJ) and writing the right side as z&(z, 2, 1), we find that 
the integral 



1 r 1 

$(2, 2, 1) = - / 1, l)dw = -Li 2 (z) 

2 Jo z 



is z^ 1 times the classical dilogarithm function. The identity (J7|) leads by induction to 
similar formulas for $(2;, 1, u) and Q(z, 2, w) when u — 2,3, ... . 

The following three examples give evaluations of <9$/<9s that we will need. 

Example 2.1. Differentiating the relation (see, for example, ^Hj) 

(16) COO = (1 - «*) 



2\ 

at s — —1, and using © and the values £(— 1) = —1/12 and 

(17) C'(-l) = ^-lnA 

where 

l^^-'-n" 

A = lim 



is the Glaisher-Kinkelin constant (see section 2.15]), we get 
18 — -1,-1,1 = In 



9s v ' ' > 2 1 /3 e V4- 

Example 2.2. Differentiating (fT7)|) and the functional equation of the zeta function [23J 
p. 27] 

C(a) = 2(27r)- 1 r(l- S )8in T C(l-a) 
at s = —2, and using Q and the value £(— 2) = 0, we deduce that 

(19) ^(12 n- 7C(3) 

(19) ^ (_1 ' _2 ' 1) -1^' 

where ((3) is Apery's constant. 

Example 2.3. Differentiating (jlj) and the functional equation of the Dirichlet beta func- 
tion [24, p. 27], 

s-l 



/?(,) = (-) r(l - 8) COS y /3(1 - S ) 
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at s = — 1, and using the value /3(— 1) = 0, we obtain the evaluation 

d® ( 1\ G 

where 

(21) G = 0(2) = i$ (-1, 2, = 2, 1)) 
is Catalan's constant. 

Theorem 2.1. For a// complex s, and complex z with §l(z) < 1/2, 

oo s \ n n s \ 

(22) (1 -*)<&(*,«,«) (^J )(«+*)"• 

n=0 ^ ' fc=0 ^ ' 

Proof. The right side of ()22j) defines a holomorphic function of z on the half-plane 3l(z) < 
1/2, because the inequality implies that |tEt| < 1- We prove ()22)1 when \z\ < 1/2, and the 
result follows by analytic continuation. 

Reversing the order of summation (for justification, see P), the right side of (j2^j) is 
equal to 

E E (I) = E (^7(1 - *)<-*>' - (i - *>•<*, ••-). 

fe=0 V 7 n=k \ / V / fc=0 V y 

and the theorem follows. □ 

Remark. Theorem 12. II excludes the Hurwitz zeta function ((s,u) = $(l,s,w). However, 
the two cases u = 1 and u = 1/2 can be obtained via the relations (J2J), ©, (QUI); and 
C (s, 1/2) = (2 s — l)C(s; !)• I n the case u — 1, formula (}22|) with z = —1 gives a globally 
convergent series for the alternating zeta function (J3J). Then (116)1 yields a formula for the 
Riemann zeta function, valid for all s ^ 1, 

CM = E^Et-D'U (* + «-. 

n=0 fc=0 ^ ' 

which was conjectured by Knopp, proved by Hasse |T2J p. 464], and rediscovered by the 
second author |T§] . 

Theorem 2.2 (Hasse). For z = 1, the series (0) extends to a function 

(23) *(!,,,„)= £ + 

n=0 fc=0 V 7 

defined for all complex s / 1, mi/i a simple pole at s = 1. Moreover, the series in (J23J) 
converges uniformly on compact sets in the s-plane to an entire function. 



Proof. We give a sketch; for details, see 

First assume that 3?(s) > 1. Setting z = 1 in (flUjl. we rewrite it as 



$(1, s, «) = ^ / -^e"-f " 2 rft = £ I 1 * J e-"V- 2 dt. 

r(s) y i-e * r(s) y ^ n + i 
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Interchanging the integral and sum, we expand the binomial and get 

Now (with k in place of n, and s — 1 in place of s) and the relation T(s) = (s— l)r(s— 1) 
yield (J2HJ) when 9?(s) > 1. 

For any s G C, the inner sum in (J2H|) with n > evidently is equal to 

(24) V(-l) fe (f) (n + kf- s = (s - 1)„ / • • ■ [ (u + x x + • ■ ■ a^) 1 -*-"^! • ■ ■ dx n , 
k=o Jo 



where the Pochhammer symbol (a) n stands for the product a(a+l) ■ ■ ■ (a+n—1). It follows 
(see [12] or 20 ) that the series in (J2"3*J| converges as required. By analytic continuation, 
this proves the theorem. □ 

Example 2.4. Note that when s is zero or a negative integer, the factor (s — l) n in (|24p 
vanishes for n > m := 1 — s, so that 

(25) E("l) fc (?) (« + fc) m = (n > m = 0, 1, 2, . . . ). 

In this case the series in ()23j) becomes a finite sum, in fact, a polynomial with ratio- 
nal coefficients. Defining the mth Bernoulli polynomial by B (x) = 1 and B m (x) = 
— m<E> (1, 1 — m, x) for m = 1,2,..., we thus have the explicit formula 



n=0 fc=0 V 7 



v3 m (x) 
In particular, 

(26) fii(z) =-$(1,0, a) =*-\ 

Example 2.5. Taking z — — 1 and s = — m in (f2*2*j) . we use (f23j) to get the explicit formula 
for the mth Euler polynomial (JT4~j) 



^) = E^E(- 1 )*(j( aj+ *) T 

n=0 fc=0 v ' 



The next two results relate the analytic continuation of $(z, s,u) for z G C - [l,oo) 
( Lemmas 12.11 and 12. 2|) to that for z = 1 (Theorem 12. 2|) . (Note that $(z, s,«) zs noi 
continuous in z at z = 1: see JT|, (fTSJl and [3 p. 30, equation (12)].) 



Corollary 2.1. For all complex s, 

v $( -z,s,u 



tfe = s$(l, s + 1,m). 



1 + 2 

Proof. Replace z with — 2 in (}2~2*|) . multiply by (1 + z)~ 2 , and integrate from 2; = to 00. 
Comparing the result with (}2*3*)) gives the corollary. □ 
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Lemma 2.3. If either z G C — { (— oo, — 1] |J [l,oo) } and s G C, or z = ±1 and 
s G C - {1}, i/iera 



(27) ${z,s,u)=T 



Proof. If \z\ < 1 (respectively, 2; = ±1 and dt(s) > 1), formula (}2"Tj) is proved by splitting 
the series (0) into two sums, one over even numbers and the other over odd numbers, and 
factoring out 2~ s . The result then follows by analytic continuation in z (respectively, in 
s). □ 

As an application, (|2"7|) with z = ±1 and s = 0, — 1, —2, . . . implies the classical relations 
among the Bernoulli and Euler polynomials |2"H equations 19:5:6 and 20:3:3]. 

Example 2.6. Differentiating (J2"?j) with respect to s at s = 0, and using (J26j) and the 



formula |p(l,0, u) = In (see p. 26]), we get 



(28) ^(-1, (),«)= In; Uj 



<9s v ' ' 7 r(^) v / 2' 
3. Double Integrals 

Using the integral analytic continuation of Lerch's function in Lemmas 12.11 and 12 .2\ we 
evaluate certain types of double integrals. 

Theorem 3.1. If either z G C — [1, 00) and 5R(s) > —2, or z = 1 and 3?(s) > —1, then 

Jo Jo 1 ~xy z u-v 

"1 /■! f^/Nu-l 



(30) / / (xy)" ^_ \ nxy ydxdy = F(s + 2)$(z, 5 + 2, iz). 

7o 7o 



Proof. The integrals define holomorphic functions of z and/or s under the conditions 
stated. Making the change of variables x = X/Y, y = Y in integral it is equal to 

lnX) s / y*" u - 1 dydX= / (-lnX) s dX 



1 — X2 7jf u — v Jo 1 — 

Substituting t = — lnX and using (flQj) (with s + 1 in place of s), we obtain equality (|2l 
when 5J(s) > and either |z| < 1 or 2 = 1; the result then follows by analytic continuation 
in z or s. To prove (jHD|), let u — > w in (f2T)j) . and use the identity (JOJ). □ 

Example 3.1. Setting z = -i, it = 1 and s = in pljl . we use (0), JUJ), fpT]) . and 
Euler's formula £(2) = 7r 2 /6 to get 

1 rl 1 vr 2 i 
-.dxdy = G 



./o 



xyi 48 
The next three examples use the substitution x — > x 2 , 1/ — > i/ 2 . 
Example 3.2. Setting z = —1, it = 1, i> = 1/2 and s = 1 in (J29|) . we get 

—dxdy — G — — . 



./o 



1 + xV " 48 
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Example 3.3. Setting z = 1, u = 1, v = 1/2 and s = 1 in (j29"|) gives 

Jo Jo 



-xlnxy , , 7r 2 
-dxdy 



\-x 2 y 2 """ 12" 

Corollary 3.1. // either z e C — [1, oo) and ri > —1, or z = 1 and n > 0, t/ien 

'- 1 (-hixy) n , , (n+l)!Li n+2 (z) 
dxdy 



o jo 



1 — xyz z 



In particular, for z E C — [1, oo) 



dxdy - ' 



o 7o (! -^)lnxy * z 
Proof. Let it = 1 in (j3~Uj) and set s = n. Then use © and (fT3j) . □ 

Example 3.4. Take z = 1/2. Using Euler's and Landen's values of the di- and triloga- 
rithm at 1/2 (see pp. 1-2]), we see that (compare Example 13.13)) 

1 rl -1 f 1 f 1 1 7T 2 In 2 2 

dxdy = In 2, / / dxdy = — 



o Jo (2-xy)\nxy J J 2-xy 12 2 

1 Z" 1 - \nxy , , 7C(3) 7r 2 ln2 In 3 2 
„ X = + — ■ 

Example 3.5. Take z = tp" 1 , <p~ 2 , — and — t^ -1 , where y> = (1 + v / 5)/2 is the golden 
ratio. Using the relation 1 + ip = ip 2 to simplify ln(l — z), and using Landen's values for 
Li 2 (z) and Li 3 (^~ 2 ) (see equations (1.20), (1.21), and (6.13)]), we obtain 

f 1 f 1 —1 f 1 f 1 1 7T 2 

-dxdy = 2 ln<£>, / / dxdy = In 2 */?, 



Jo {(p-xy)lnxy ' r ' J Q J <p - xy 10 

1 rl i /■! X It 2 

dxdy = In </?, / / — oxo 1 ?/ = In 2 ip, 



o Jo (<p 2 -xy)\nxy ' r " 7 7 <p 2 - xy ' 15 

dxdy = — , f f dxdy = h r 



o Jo (l + <pxy)lnxy ' <p J J l + cpxy ' lOy? 

/ / 7 — : Yi dxdy = \ncp, [ [ — ^ drrdw = ^- - 

/o </o (^ + xy)lnxy J J tp + xy 15 2 

1 Z" 1 -lnxy , , 8C(3) 4vr 2 ln^ 41n 3 <^ 
o Jo V 9 - xy 5 15 3 

Example 3.6. Using Ramanujan's values jHl Part IV, p. 324] for Li 2 (— 1/8) +Li 2 (l/9), 
Li 2 (— 1/2) + |Li 2 (l/9), and Li 2 (l/3) — ^Li 2 (l/9), and substituting x, y — > x 2 ,y 2 in our 
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double integral for the last Li 2 (l/9), we get, respectively, 

/ / (8 + xy)(9-xy) dxdy= 2 ln 8' 
f f 1 $2-7xy = vr 2 _ 

f 1 f 1 9 + xy , , vr 2 In 2 3 

Example 3.7. In Bailey, Borwein, and Plouffe's [Tj dilogarithm "ladder" 6Li 2 (l/2) — 
6Li 2 (l/4) - 2Li 2 (l/8) + Li 2 (l/64) = C(2), we replace Li 2 (l/2 fc ) with a double integral in 
which we substitute x, y — > x k , y k , for k = 1, 2, 3, and 6, obtaining 

f 1 f 1 (4 + 2xy + x V) (4 - 8xy + x 2 y 2 ) , , vr 2 

I I ^¥ dxdy = n 

Corollary 3.2. We have 

(31) I' C ^ lnXy ^ dxdy = T(s + 2)((s + 2) (&(s) > -1), 
Jo Jo 1 ~ X V 

(32) T ("^^' dxdy = r( S + 2)C(5 + 2) > -2), 
Jo Jo 1 + #2/ 



(33) r T ( l^' dsdy = r( S + 2)/?(s + 2) > -2). 

Jo Jo 1 + ^ y 

Proof. In ()30j) take z = 1, it = 1; z = —1, it = 1; and z = —1, u = 1/2; then use relations 
(J2J), © an d (jU), respectively. After using (J1J), substitute x, y — > x 2 ,y 2 . □ 

For integer s > 0, formula (J3T|) is due to Hadjicostas 

Example 3.8. Taking s = and 1 in (|31|). we recover Beukers's integrals [3] 



(34) 



/ / dxdy = ((2), [ [ — dxdy = 2((3). 

Jo Jo 1 - X V Jo Jo 1 - X V 



Example 3.9. Taking s = —1, 0, and 1 in (|33|) . and using (J2TJ) and the values = 7r/4 
and /3(3) = vr 3 /32 (see [21 p. 29]), we get 



— 1 . . 7T / / 1 



'0 JO 



dxdy = -, / / — — —dxdy = G, 



[1 + x 2 y 2 ) \ia xy 4' y 7 1 + x 2 y 2 

1 — hixy 7T 3 

o Jo l + ^V 16 

Example 3.10. Let p be a non-trivial zero of the zeta function (J2J). Differentiate (|32j) at 
s = p — 2 and use (fTo^) to get (compare Example 14. 3|) 

' ' 1 y; -ln(-lna;y)da;dy= (1 - 2 3 ^)r(p)C'(p)- 



/o ./o 



1 + xy 
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Corollary 3.3. If either z G C - ((— oo, — 1] U[l, oo)) and dt(s) > -2, or z ± 1 and 
$t(s) > —1, then 



1 r> '-\nxy) s ^Xs+2(z) 
Y^dxdy = T(s + 2) 



Jo 



1 — x^y^z* z 
where Xs(z) ^ s Legendre's chi function 

Proof In flHUj) . take u — 1/2 and substitute x, ty, z — > x 2 , y 2 , z 2 . Then use (J5J). □ 

Example 3.11. Set s = 0, and take z = tan | = \/2 — 1 and 2 = <£>~ 3 = — 2. Using 
Landen's values [TBI section 1.8.2] for x 2 (tan|) and x 2 (<^~ 3 ), we obtain, respectively, 



1 r i 



^0 



1 



1 — x 2 y 2 tan z § 



—dxdy 



7T 



In 2 tan ? 



1 a-i 



/o </o V - 
Corollary 3.4. For z G C — [1, oo), 



x 2 ?/ 2 



dxdy 



16 tan f 4tanf ' 
7T 2 3 In 2 p 



24(^ 3 



4yj 3 



(35) 



(36) 



1 /■! 



1 — £2/2) In a;?/ 
-(xy) 



dxdy 



./o 

1 /•! („„.\u-l 



U — V 

dxdy = $(2, 1, w 



«9$ 9$ 

^(^0,,)-^(^0,n) 



./o (l-a^lnacy 

Proo/. Substitute (s + l) _1 r(s + 2) for T(s 4- 1) in and let s -> -1+. Since 2 G 
C — [1, 00), in view of (fT2"|) the result is (f3"H|) . Letting t> — > u in (|33j) . and using identity 
©, we arrive at flUJ). □ 

Example 3.12. Set z = 0. Using (JTJ), we get 

•1 rl /„..\u-l 



^0 



dxdy 



J- in" 



it — f w 



(xj/T^ 1 

— axay = — . 

o — In xy it 



— In xy 

Example 3.13. Set z — 1/2, u — 2, and v — 1 in 1)35)1 . Using the derivative of (0) with 
respect to s at s = 0, we get (compare Example 13 ,4jl 



1 rl 



-X 



Jo (2-xy)lnxy 



dxdy = In a, 



where 
(37) 



a 



l\/2v / 3^=l 1/2 2 1 / 4 3 1 / 8 --- 



exp 



19$ 
2 ~ds 



,0,1 



is one of Somos's quadratic recurrence constants ^H] (see p. 446], which uses the 
notation 7 instead of a; see also [TOJ P- 348]). 

Corollary 3.5. For u > and t> > 0, 



(38) 
(39) 



1 rl 



JO 

7" 

^0 



'1 + xy) Inxy 
-(xy) 11 - 1 , , 1 

; — axctu = - 

1 + xy) mxy 2 



1 r(-) r(^±i 

dxdty = In 



u — V 



r(|) r(2±i) 



u + l 
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where ip(u) = 4-lnT(u) = T'(u)/T(u) is the digamma function. 

Proof. To prove (JHEJ), set z = —1 in (JH5J) and use (J2BJ)- To prove (JHIIJ), let v — > u in 

(JSHJ)- □ 

Example 3.14. Setting u = 1 in (|39J) and w = 2 and t> = 1 in (J3*%|) gives integrals for 
In 2 and ln(7r/2), respectively. Their sum and difference are integrals for ln-zr and the 
"alternating Euler constant" ln(4/7r) (compare Example 14.41 and see and |23j): 

1 rl 1 + x , , f 1 f 1 l-x , , , 4 



dxdy = Inn, / / — — — -dxdy = In — . 



/o Jo (! +a;y)(-lna;y) ' J J (1 + xy)(-\nxy) ' 7T 

Example 3.15. Set it = 1 and v = 1/2 in (J3~SJ). and substitute x — > x 2 , y y 2 . Using 
the reflection formula r(t)T(l — t) = ir/ sin7rt with t = 1/4, we get 

1 rl -x , . ^Phx 



-dxdy = In ■ 



/ J (l + x 2 y 2 )\nxy^» '"r(3/4) 2 ' 
Corollary 3.6. For u > and v > 0, 

(40) r /' - *'<") - 



o Jo 



1 — xy u — v 

(41) f f { ^-dxdy = ^\u). 

Jo Jo 1 ~ • L y 

Proof. To prove (JjDJ) , set z = 1 in (J29|) . let s — > 0, and use the formula (23 p. 271] 

(42) lim s, u) —J = 



*-»i V s- 1, 

To prove (JUJ), let v -> u in (JUJ). □ 

Example 3.16. Take w = 1/3 and v = 2/3 in (HOJ) . Since ^(2/3) - ^(V 3 ) = (see 
equation 44:7:1]), the change of variables x — * x 3 and y — > yields 



o 7 1-xV 3^3" 

Example 3.17. Taking u = 1/2 in (JUJ) and using the value ^'(V 2 ) = 7r 2 /2 (see [2U p. 
427]), we substitute x — ► x 2 , y — ► ?/ 2 and obtain 



1 .1 1 ^2 



1 — X 2 ?/ 2 



dxdy 



>o Jo 

Corollary 3.7. IfU(s) > -2, then 

(43) / I x u - 1 y v - 1 (-\nxy) s dxdy = r(s + l)-' 



o 



1J — If 

„u-i„.v-i,_ i nxv y dxdy = r ( s + !) 

W — V 

/•I 



(44) I [ (xy^i- In xy) s dxdy = r(s + 2)u 

Jo Jo 

Proof. Set z = in (J2H and (JHJ), and use ©■ □ 
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Example 3.18. Taking s = —3/2, u = 2 and v = 1 in (jj^) gives 

Example 3.19. Taking u = 1 and s = —3/2 or s = —5/4 in (|44|) gives 
1,1 1 ,1,1 1 




dxdy = y/n, I I - — -^j-dxdy = T (3 / A) . 



/o Jo (-\nxyf 2 " V ' Jo Jo (-Inxy)^ 
Corollary 3.8. I/z G C - [1, oo) and 3fc(a) > -2, tfien 

(45) /' / 'fS'-'* 

= - u)$(;z,s + l,v) + $(z,s,<u) + (it - l)$(z,s + l,u) - $(z,s,u)], 

( 46 ) f [ 77^ ^(-\nxy) s dxdy = T(s + 2)1(1 -u)$(z,s + 2, u) + <S>(z,s + l,u)]. 

Jo Jo (1 - xyz) 2 

Proof. Differentiate ()29|) with respect to z and apply identity (JSJ). Replacing w with w — 1, 
and v with t> — 1, we arrive at (JJ3J). To prove (j4*Uj) . let w — > u in fj43j) and use identity 
©• □ 

Example 3.20. Set z = —1, u = 1/2, i> = 1, and s = 2 in (|4~5j) . Substituting x — > x 2 , 
y -> y 2 , and using ©, ©, (ED, and the values C*(2) = vr 2 /12 and /3(3) = tt 3 /32, we get 

' xln xy 7r 7r° 

Example 3.21. Set z = —1, it = 1 and s = 4 in (|4l)j) to get 

dxdy = —— C(5). 





Example 3.22. Setting z = —1, it = 1/2 and s = — 1 in (JIEJ), and substituting x — > x 2 
y -> y 2 , we get 

So So (l + x 2 y 2 ) 2 \nxy dxdy = "IT" 
Corollary 3.9. I/z G C — [l,oo) ; £/ten 



1 ,1 _ x U—lyV—l 

dxdy 




o ./o 



(1 — xyz) 2 \iaxy 
1 

— f i - >')■ , ,, , , . , , , , ,, , 

OS OS OS OS 



(9$ (9$ 9$ 9$ 

(1 - v)—{z, 0, v) + — (z, -!,«)-(!- u )— (z, 0, u) - j-(z, -1, u) 



u — V 

Proof. Substitute (s + l) _1 r(s + 2) for T{s + 1) in gSJ) and let s -> -1 + . □ 
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Example 3.23. Setting z = -1, u = 2 and v = 1, we use PI, £EEJ), and (0) to get 
(compare [221 equation (9)]) 



1 r i 



—x 



'o (l + xy) 2 \nxy 
where A is the Glaisher-Kinkelin constant. 



dxdy = In 



A 6 



2V8 



Example 3.24. Set 2 = —1, w = 2 and i> = 1/2. Using (|2lljl and substituting x — > x 2 
and y ^ y 2 , we get 



1 /•! 



/ / I ., — dxdy = — . 
o Jq (1 + x z y z ) z mxy n 

For the proof of the final theorem in this section, we require a lemma. 

Lemma 3.1. If H n r = 1 + ^ + • • • + ^ is the nth harmonic number of order r and 
H n = H n ^i, and if \z\ < 1, then 



(47) 



00 



+ 2 V ^z n = 3Li 3 (z) - Li 2 (z) ln(l - z). 



n=l 



77 



n=l 



Proof. We have 



00 n—l 



00 00 



:v :v z k+j 



H „2 z Li3 W SS n2A , E E n 2 fc EE(h /)"/, 

n=2 fc=l fe=l n=fc+l fe=l j=l v J/ 



n=l 



E7E/ 7/ «' !+j - 1 <i«* = £ 

fc=l ft ,=1 ^0 1 J o fc= i 



1 p 1 



dudt 



1 /"'-ln(l-M 



^ Jo 



dudt 
l-u 2 



1 fln^l-t) 



and in a similar way 



E^V-LisM 

n=l 



Li 2 (*) 



1- t 



Adding twice the first identity to the second, we obtain 



V^" + 2V^™-3Li 3 (^ 



Li 2 (t) + ln 2 (l-t) 



1-t 



7j 



dt = -Li 2 (z) ln(l - z) 



n=l n=l 

the evaluation of the integral being easily checked by differentiation. 
Theorem 3.2. For z G C — [1, 00), i/ie following analytic continuations hold: 
1 H ln(l-xz) 



□ 



(4 



(49) 



J 



7 

./ 



1 — xy 
'1 



-dxdy 



lnzln(l -z) + Li 2 (l - 2) 



ln(l — xyz) 
1 — xy 



ln(l-z)+Li 3 (z)-Li 3 (l-z) + C(3), 



dxdy = Li 2 (l - 2) ln(l - z) - Li 3 (z) - 2Li 3 (l - z) + 2((3). 
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Proof. By analytic continuation, it suffices to prove (|48|) and (|49j) when \z\ < 1. 

Expanding — ln(l — xz) in a series, and using (|4Hj) and the identity ^(n+l) — ^(1) = J/r, 
(see j23 equation 44:5:5]), we see that 



[ [ — -7— —dxdy = Y] — [ [ dxdy = V — ' 

Jo Jo 1 - ^ f-f « Jo Jo 1 - x v r-f n 



71=1 »U n=1 



Ramanujan's summation of the last series p. 251, Entry 9(i) and p. 259, equation 
(12.1)] completes the proof of (EEHl) . 

Similarly, using (}4~Tj) and the identity if)'(n + l) = \n 2 — H n 2 (see j2U equation 44:12:6]), 
we find that 



Jo Jo i-xv 6 1 ti n 



Now apply (@ZD with J2n=i H n n~ 2 z n replaced by the right side of ()48j) . Simplifying with 
Euler's identity Li 2 (z) + Li 2 (l — z) = |7r 2 — lnzln(l — z) (see [HI p. 1]), we arrive at 
(|49j). □ 

Example 3.25. First take 2 = 1/2; write In (l — |) as ln(2 — x) — In 2 in (jlSjl. rewrite 
f!49jl similarly, and apply the first equation in ()34j) . Then take z = —1. Using the values of 
Li 2 (l/2) and Li 3 (l/2) (see Example 13 .5|) . and the inversion formulas for Li 2 (y) + Li 2 (l/y) 
and Li 3 (y) — Li 3 (l/y) (see equations (1.10) and (6.7)]) with y = 2, we obtain 

^ dxdy = I ( — ^ — — — —dxdy= — C(3), 



Jo Jo 1 - xy Jo Jo 1-xy 8 



l-xy J J 1 - xy 4 



Example 3.26. Letting z — > 1 , we deduce that 

-1 rl ™N rl rl 



I f -^ l - X h xd y = 2/ / -■°< 1 -^^ y = 2 C (3). 
Jo Jo 1 - ^2/ Jo Jo 1 - xy 

4. A Generalization of Hadjicostas's Formula 

Hadjicostas ^1] asked for a generalization of the double integral formulas in JUj and 
for the Riemann and alternating zeta functions ( Examples 14.11 and 14. 2|) . To fill the 
bill, we multiply the integrand of (j3*Uj) by 1 — x, and obtain a formula valid on a larger 
half-plane of the complex s-plane. 

Theorem 4.1. If either z 6 C — [1, 00) and 9R(s) > —3, or z = 1 and 9ft(s) > —2, then 

(50) / I l ~ X ( xyY' 1 (- In xyY dxdy 
Jo Jo l-xyz 

[1-Z)$(Z 1 8 + 1,U)-U- S - 1 



r (s + 2) 



®(z,s + 2,u) + 



z(s + l) 
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Proof. The integral in (}50|) defines a function which is holomorphic in s, when 3?(s) > —3 
if z e C - [1, oo), and when 3?(s) > -2 if z = 1. We prove (J501) with 3?(s) > and the 
result then follows by analytic continuation. 

Replace u with u + 1 in (j2HJ) and set v = u. Subtracting the result from (|3T)J) gives 



(51) 



1 rl 



J 



— In xy) s dxdy 



1 — rq/z 

= r(s + 2)$(z, s + 2,u) + r(s + 1) s + 1, u + 1) - s + 1, u)] . 
Using identity (J7J), we obtain ([50]). and the theorem follows. □ 
Example 4.1. The case z = w = 1 is Hadjicostas's formula [S], [TUj (see also 



(52) 



1 r l 



JO 



X 



1 — xy 



hixy) s dxdy = T(s + 2) 



C(s + 2) 



8 + 1 



m*) > -2). 



Example 4.2. Taking z = — 1 and u = 1 gives the analogous formula 



1 rl 



1 — X 

o Jo l + zjj 



\iaxy) s dxdy = F(s + 2) 



C( S + 2) 



l-2C*(s + l) 

8 + 1 



(3?(s) > -3). 



Example 4.3. Let p be a non-trivial zero of ({s). Differentiating (|52jl with respect to s 
at s = p — 2, we get (compare Example 13.10)1 

-1 /•! 



o ./o 



1 — xy 

Corollary 4.1. For u > 0, 

-1 /•! 



(— \nxy) p 2 ln(— In xy)dxdy = T(p) 



C(p) + 



1 



<hp) 



(p 



p 



o 



1 — X 



(xy) u dxdy = \nu — ifj(u). 



□ 



o (1 -xy)(-\nxy) 
Proof. Set z = 1 in PH), let s -> -1+, and use IfEZj) . 
Example 4.4. Setting u — 1 and using the value |2*H p. 427] 
(53) = -7, 

+ - — Inn) is Euler's constant we recover the formula in 



where 7 = lim^oo (l + | + 
and 



1 /•! 



jo 



1 — x 



[1 -xy)(-\nxy) 



dxdy = 7. 



Corollary 4.2. For z G C — [1, 00), 



1 r l 



./0 



1 — X 



1 — xyz) In xy 



[xy) u dxdy 



= —(z,0,u) + - 
os z 

Proof. Let s -> -2+ in flHU) and use (THZJ) . 



9$ 

(2; — 1)— — (z, — 1, w) — w lnw 
as 



+ 



□ 
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Example 4.5. Take u = 1/2 and z = — 1. Substituting x — > x 2 , y — > y 2 and using (|2 
and pOjl . we get 

1/11 1-x 2 , , , r(l/4) 2G 1 

axay = In 



/ 7o (1 + xV) In 2 2r (3/4) tt 2' 

Example 4.6. Replacing $ with the series (pQ), we let z — > and deduce that 

/ / ^^(xy)"- 1 ^?/ = (u + 1) In fl + -1. 
Jo Jo In xy V «/ 

Corollary 4.3. // 9?(s) > -3, t/*en 

/ / (l-a;)(a;y)"- 1 (-lnx?/) s rfxrf?/ = r(s + l) [(s + 1-u)m- s - 2 + (m + 1)- 5 - 1 ] . 
Jo Jo 

Proof. Set z = in flSU) and use ©. □ 
Example 4.7. Taking s = —5/2 and w = 1 gives 

rr =3 r (8 ^- io) - 

5. Infinite Products 

Using corollaries of Theorems 12.11 and 12. 2^ together with logarithmic series for the 
digamma and Euler beta functions (Theorems 15.11 and 15.2)1 . we derive infinite products 
for many constants. 

Corollary 5.1. For m — 0, 1, 2, . . . and complex z with dt(z) < 1/2, 

v 7 n=0 v 7 fc=0 V 7 

Proof. Differentiate ()22)1 with respect to s at s = and multiply by (1 — z)~ x . Then apply 
(JHJ) m times. □ 

Example 5.1. Setting m = 0, z = — 1 and it = 1 and multiplying by 2 gives 

n=0 k=0 v 7 

Using (J2HJ), we recover the product from 

tt /2\ 1/2 / 2 2 \ 1/4 /2 3 -4\ 1/8 / 2 4 -4 4 xl/16 



2 v 1 / v 1 " 3 / v 1 " 33 / V 1 - 36 " 5 , 

Example 5.2. Setting m = 1, z = — 1 and u = 1 and multiplying by 4, we get 

^(-i.-m)=E^D-i)-@m» + i). 

n=0 fc=0 v 7 
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Using (fl8|) . we obtain the product 

A i2 / 2 \ 2/2 f 2 2 \ 3/4 /2 3 -4\ 4/8 / 2 4 -4 4 \ 5/16 

which converges faster than the products involving A in Examples 15.71 and 15.111 
Example 5.3. Setting m = 2, z = — 1 and u — 1, we get 

n=0 fc=0 v 7 



Then ([IS}) gives 



3 7C(3)/4tt 



2 



1/8 / 2 \ 3/16 /o3 „\ 6/32 / 2 4 • 4 4 ^ 10/64 



1-3/ V 1 • 3 3 / Vl-3 6 -5 



Example 5.4. Take m = 0, z = —1 and w = 1/2. Writing In (k + |) as \n(2k + 1) — In 2 
and using (|23|) . we see that 



ds \ • "2 2 



oo _ n ✓ \ 

E^E(- 1 )' +1 u)M2* + i). 

»i— n i — n V / 



n=0 fc=0 

Using (|2*8j) . we get the product 

T(l/4) /3\ 1/4 / 3 2 \ 1/8 /3 3 -7\ 1/16 / 3 4 ■ 7 4 ^ 1/32 



2r(3/4) V 1 / V 1 " 5 / V 1 " 53 / V 1 " 56 " 9 
Example 5.5. Setting m = 1, z = —1 and ii = 1/2, we get 

v 7 n=0 fc=0 v 7 



Use (J2U)l and exponentiate to obtain 

3 G/tt _ 



e 



3 \l/8 / 32 \2/16 / 3 3. 7 \ 3/32 , 3 4 . 7 4 \ 4/64 



\) V 1 ' 5 / v 1 " 53 / v 1 " 56 " 9 

Corollary 5.2. For a// complex s / 1, 

4a , s , u)+(3 _ 1) ^ (M , uH |^|: ( _^( : )M^. 

Proof. Multiply ()23|) by s — 1 and differentiate with respect to s. □ 

Example 5.6. Set u — 1 and s = and multiply by 2. The relation (J2J) and the values 
((0) = -1/2 and C'(0) = -\ In 2tt give 

oo n , x 

In 2tt - 1 = £ — 5>l) fe+1 (")(* + 1) M* + !)• 

n=0 fc=0 ^ 7 

Exponentiating, we obtain the product 

2tt /2 2 \ 2/2 / 2 4 \ 2/3 / / 2 6 -4 4 \ 2/4 / 2 8 -4 16 \ 2/5 



1) V 1 • 33 / V 1 • 39 / V 1 • 318 • 55 



18 JESUS GUILLERA AND JONATHAN SONDOW 

Example 5.7. Set s = — 1 and u = 1 and multiply by 1/2. Using relations (J2J) and fTTj) . 
and the value £(— 1) = —1/12, we get the formula 

1 °° i n / \ 

(«) "-r^ Ei- 1 )" 1 Q <* + d 2 >»<* + 1). 

n=0 fc=0 V 7 

which yields 

A /2 4 \ 1 ^ 4 / 2 8 \ /2 12 -4 16 ^ 1 ^ 8 



eV8 V 1 / V 1 • 3 9 J V 1 • 3 27 
Theorem 5.1. For u > 0, 

oo 1 n r \ 

(55) V(w) = E ^TT J ln (« + *)■ 

n=0 fe=0 ^ 7 

Proof. Multiply ([2*3*]) by s — 1 and compute the derivative with respect to s at s — 1, using 

□ 



Example 5.8. Take u = 1, multiply by —1, and exponentiate. Using ([53)1 . we recover 
Ser's product [T7j (rediscovered in [21] and [215] ) 



e 7 



1/2 / 2 \ V3 / 2 3 . 4 \ 1/4 / 2 4 . 4 4 \ 1/5 



1-3/ V 1 " 33 / V 1 " 36 " 5 
Compare the remarkably similar product for e in Example 15.121 

Example 5.9. Take first u — 1/4 and then u = 3/4. Using the relation if) (3/ 4)— if) (1/4) 
n (see [2H p. 427]), we obtain the series 



El A, n fc ^ n ^ln 4 ^ + 3 
n + 1 ^ UJ n 4A; + 1 



7T 

n=0 ' fc=0 

which gives the product 

'3 \ V3 • 5\ V3 ■ 5 2 ■ 11\ V3 ■ 5 3 ■ ll 3 ■ 13\ V3 ■ 5 4 • IP • 13 4 ■ 19 



1) \l-7j V 1 ■ 7 2 ■ 9 J V 1 ■ 7 3 • 9 3 ■ 15 J V 1 ■ 7 4 ■ 9 6 • 15 4 ■ 17 
For the next result, recall the formulas for the Euler beta function [23 p. 254], 

B(u,v)= / x (1 - x) v L dx 



T(u + v) 
Theorem 5.2. For j — 1,2, ... , 

(se) b(«, = e — n E(- 1 ) fc+1 J ln (« + k). 

n=j J k=0 ^ 7 

Proof. We first establish the formula 

(57) £*J^^ dx = ± { - ir ^ Hu + k) 

JU fc=0 v 7 

for n — 1, 2, . . . . Fix n and define /(«) to be the difference between the integral and 
the sum in ([57]) . Then /'(w) = 0. (To see this, differentiate f(u) under the integral sign, 
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expand (1 — x) n , and evaluate the integral.) Thus f(u) = C is a constant. By (}25|) with 
m = 0, we may replace ln(w + k) with In (l + -) in /(w). It follows that f{u) — > as 
m — > oo. Hence C = and (|H7j) follows. Now to prove (JHEJ), multiply (|H7j) by (n — j + 
and sum from n = j to oo. □ 

Example 5.10. Set u = 1 in (|H7j) . multiply the equation by (— l) n+1 , and sum from n = 1 
to oo. Using Example 13.131 (let x = X/Y, y = Y and integrate with respect to Y), we 
obtain 

where a is Somos's constant (jSZj). Exponentiation gives a product for o with rational 
factors: 

_ 2 1 ■ 3 2 3 ■ 4 1 • 3 6 ■ 5 
(58) ° ~\' ~¥~ ' 1^33 ' "2471^"' ■ 

Since the partial products are alternately above and below a, and converge slowly, the 
product is a good candidate for acceleration by Euler's transformation [T5], [23! • Applying 
it and simplifying (as in [231 section 3]), we arrive at the product 

2 \l/2 / 3 \V4 / 4 \ 1 /8 / 5 \1/16 



(59) o- ,,,,,, 

l iy v 2 7 V 3 7 V 4 , 

which converges faster than the products for a in both (|3T|) and (J3HJ) - (Note that 
also follows by writing a as <t 2 /o" and substituting the product in (|3T)l .) 

Example 5.11. Set it = 1 and j = 4 in 1)56)1 . and replace n with n + 4. Multiplying the 
equation by 1/2, so that the left side becomes 5(l,4)/2 = 1/8, we add the result to (jSljl . 
Replacing with k — 1, and n with n — 1, we obtain 

00 1 r n y 1 v n+4 / i Q\ 

n=i Lfe=i v 7 fe=i v 7 

This yields a product for the Glaisher-Kinkelin constant: 

.4 



2 4 . 4 4 \ ^ 2 /2 9 ■ 4 10 ■ 6\ 1//4 / 2 14 ■ 4 20 ■ 6 6 x 



1 ■ 3 6 ■ 5 y Vi-3 10 -5 5 y V 1 " 324 - 515 - 7 

The remaining results use only the series for the Euler beta function in Theorem 15.21 
In particular, their construction does not involve Lerch's function. 

Example 5.12. As B(l,j) = 1/j, exponentiating (|56jl with u = 1 yields a product for 
e 1 /- 7 . When j = 1, this is the first author's product for e in [23] : 

e = 



2 \i/i / 2 2 X 1 / 2 /2 3 -4\ 1/3 / 2 4 -4 4 N 1/4 



The extension to e 1//j for j = 2, 3, . . . was first found by C. Goldschmidt and J. B. Martin 
[S] using probability theory. With j = 2, it gives 

,2 \ 1/1 / 2 3 . 4 \ 1/2 / 2 4 . 4 4 N 1/3 



e 



1-3/ Vl-3 3 / V 1 • 3 6 - 5 
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Theorem 5.3. The infinite product representation of the exponential function 

l/n 



oo In 



(60) e x = J] \[{kx + 1) W 

71=1 \k=l 

is valid for any real number x > 0. 

Proof. The formula holds when x — 0. If x > 0, then since it, unlike j, is not restricted 
to integer values in Theorem 15 .2\ we may take u = 1/x. By (I25j) with m = 0, we may 
then replace In (- + k) with \n(kx + 1) in (j5fijh Setting j = 1 and using B(l/x, 1) = x, 
we exponentiate (JHHj) and obtain (JfiOj). □ 

Example 5.13. Taking a; = 2, we get 



e 2 



1/1 / Q 2 \ 1/2 / 3 3 . y\ 1/3 / 34 . 7 4 \ 1/4 



1-5/ V 1 " 53 / v 1 " 56 - 9 . 
When x = p/q is a rational number with p > and g > 0, we can simplify ()60j) by 
replacing kx + 1 with fcp + g, again using ()25j) with m = 0. An example is 



e 2/3 



1/1 / & \ V 2 /53 . g\ 1/3 / 5 4 . g4 \ 1/4 



3 ■ 7 / V 3 ■ 7 3 / V 3 ■ 7 6 ■ 11 



In particular, with x = 1/q and q = 2, 3, ... , the roots e 1//g are represented as different 
products from those in Example 15.121 For instance, 



1/1 / q2 \ 1/2 / 3 3 . 5 \ 1/3 / g4 . 54 \ 1/4 



2-AJ \2-4 3 J V 2 " 46 " 6 
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